Strong electron-phonon coupling in the rare-earth carbide superconductor La 2 C3 
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We present the results of a crystal structure determination using neutron powder diffraction as 
well as the superconducting properties of the rare-earth sesquicarbide La2C3 (T c » 13.4 K) by means 
of specific heat and upper critical field measurements. From the detailed analysis of the specific 
heat and a comparison with ab-initio electronic structure calculations, a quantitative estimate of 
the electron-phonon coupling strength and the logarithmic average phonon frequency is made. The 
electron-phonon coupling constant is determined to \ p h ~ 1.35. The electron-phonon coupling to 
low energy phonon modes is found to be the leading mechanism for the superconductivity. Our 
results suggest that La2C3 is in the strong coupling regime, and the relevant phonon modes are La- 
related rather than C-C stretching modes. The upper critical field shows a clear enhancement with 
respect to the Werthamer-Helfand-Hohenberg prediction, consistent with strong electron-phonon 
coupling. Possible effects on the superconducting properties due to the non-centrosymmetry of the 
crystal structure are discussed. 

PACS numbers: 74.25.Bt, 74.25.Jb, 74.70.Ad, 74.62.Bf 
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I. INTRODUCTION 

The recent discovery of superconductivity in MgB^ 1 
and alkaline earth-intercalated graphites 2 as well as Li 
at high pressures^ has renewed the interest in electron- 
phonon (e-ph) coupled superconductors without strong 
electron-electron correlations. The electronic states and 
phonon modes relevant for superconductivity vary from 
system to system, but in general superconductivity in 
these compounds benefits from the light atomic mass of 
constituents e.g. of boron, carbon and lithium. When a 
specific part of the Fermi surface (FS) couples strongly 
to high frequency phonon modes, an increase of T c can 
be achieved even if the e-ph coupling averaged over the 
full FS remains moderate. This mechanism opens a new 
route to achieve high T c 's by e-ph pairing. 5 . This is often 
associated as the " metallic hydrogen superconductivity" 
scenario £ 

The rare earth sesquicarbides, i?2C3 (R — Rare 
earths), which crystallize in the bcc PU2C3 structure- 
type, have been suggested as possible candidates, in 
which such conditions are realized. Early on, T c of La2C3 
and Y2C3 was found to be ~ 11 K, and Th doping in 
Y2C3 raises T c to 17 comparable with the T c 's of 
the A15 compounds. In the PU2C3 structure, C2 dumb- 
bells are located inside the rare earth metal atom cage. 
Since the C-C bond is quite short, the phonon frequency 
for the C-C stretching phonon modes is expected to be 



very high. In fact, recent ab-initio calculations showed 
that the C-C stretching phonon frequency in Y2C3 is 
^1442 cm -1 . 8 Therefore e-ph coupling between the high 
frequency phonons and C-C antibonding states at the 
Fermi level (E-p) has been considered as an origin of the 
relatively high T c in Y2C3 and La2C3. 

Superconductivity in rare earth sesquicarbides recently 
regained attention because of the discovery of 18 K su- 
perconductivity in Y2C3 samples prepared under high 
pressure (~ 5 GPa)^^ The upper critical field (H c2 ) 
is also significantly increased and amounts to H C 2(0) 
>30 T. 9 Electronic structure calculations 8,11 for Y2C3, 
however, demonstrated that the high frequency C-C bond 
stretching phonon modes contribute less than 10% of 
the total e-ph coupling, and rather low frequency Y(-C) 
phonons must be considered to be the relevant modes 
for superconductivity For stoichiometric band filling, 
the e-ph coupling constant X p h is predicted to be 0.6. 
These results were consistent with previous H C 2 studies 
on La2C3, which reported X p h ~ 0.8 suggesting moderate 
e-ph coupling. 12 However, such a moderate e-ph coupling 
appears to be too small to generate the relatively high T c 
in the sesquicarbide superconductors^ Therefore, to rec- 
oncile with the measured T c , it is necessary to consider 
either significant e-ph coupling with the high frequency 
C-C stretching phonon modes or strong e-ph coupling 
with the low frequency rare earth related phonon modes. 
In order to shed light on this controversy, further exper- 
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imental studies on the superconducting properties of the 
R2C3 are required. 

In this paper we focus on one of the rare earth 
sesquicarbides, La2C3. We investigate its superconduct- 
ing properties together with the crystal and electronic 
structures. In the rare earth sesquicarbides, it has been 
well known that T c as well as the structural properties 
vary significantly depending on the synthesis and anneal- 
ing conditions, essentially due to C deficiency. Therefore 
in order to carry out reliable studies on such compounds, 
one needs to characterize the superconducting properties 
simultaneously with the structural properties with par- 
ticular attention to C deficiency. Unlike Y2C3, which re- 
quires high pressure and high temperature preparation, 
La2C3 can be prepared at ambient pressure conditions 
using standard arc melting technique, and larger sam- 
ple quantities can easily be synthesized. Recently, we re- 
ported that T c of La2C3 can be enhanced up to 13.4 K us- 
ing excess C in the starting composition of the materials 
combined with an adequate post-annealing i 13 i 14 i 15 i 16 i 17 

In addition, La 2 C3 can be a potential system for inves- 
tigating the effect of non-centrosymmetry in the struc- 
ture. The space group symmetry of R2C3, /43d, belongs 
to the tetrahedral crystallographic class Td lacking a cen- 
ter of symmetry. When the crystal structure has no cen- 
ter of symmetry and spin-orbit coupling is significant, 
the degenerate spin-up and spin-down bands are mixed 
and split, which can induce unexpected superconduct- 
ing properties. For example, one of the well-known non- 
centrosymmetric superconductors CePtsSi^ where Ce 
/-bands possess significant spin-orbit coupling,- 9 , shows 
unconventional superconducting properties such as high 
H C 2(0) exceeding the Pauli limits and a line node in 
the superconducting order parameter. 20 ' 21 However, this 
system also shows a heavy fermion nature, and it is not 
clear yet how far the origin of the exotic properties has 
to be attributed to the non-centrosymmetry. Thus, it is 
required to explore other non-centrosymmetric supercon- 
ductors without strong electron correlations^ In this re- 
spect, -R2C3 compounds with non-magnetic rare earths, 
can also be such model compounds to study only the 
effect of the non-centrosymmetric structure without the 
interference from the magnetism of the constituents. 23 
In case of Y2C3, the spin-orbit coupling for Y d-bands 
may be rather small, thus we can not expect a signifi- 
cant effect of non-centrosymmetry. Because of the higher 
atomic mass of La as compared to Y, spin-orbit coupling 
effects in La2C3 is expected to be more pronounced than 
in Y2C3. Note that La is placed next to Ce in the periodic 
table and accordingly is expected to generate comparable 
spin-orbit coupling. 

The paper is organized as follows; we first provide ex- 
perimental details including a brief description on the 
synthesis procedures (Sec. II). Secondly we discuss the 
crystal structure at low temperatures gained from neu- 
tron powder diffraction (NPD) investigations (Sec. Ill) 
and report the electronic structures of La2C3 obtained 
from ab-initio calculations (Sec. IV). Specific heat (Sec. 



V) and the upper critical fields (Sec. VI) on the samples 
characterized by NPD are presented. Finally we will dis- 
cuss the e-ph coupling strength of La2C3 and also the 
effects of non-centrosymmetry on the superconducting 
properties (Sec. VII). 



II. EXPERIMENTAL 

Polycrystalline samples of La2C3 were synthesized by 
arc melting the constituents on a water-cooled Cu cru- 
cible under purified Ar atmosphere. La metal chips 
(Ames Laboratory, 99.99%) were used with spectroscopic 
grade graphite chips (Deutsche Carbone, 99.99%). Be- 
fore use, the graphite chips were outgased overnight at 
950 °C under high vacuum conditions (P < 1CP 5 mbar). 
The starting materials were melted more than 6 times 
and at each time the button was turned over to ensure 
homogeneity. La2C3 is very moisture-sensitive, tending 
to decompose readily within a few minutes on exposure 
to the air. Consequently all handling of the starting ma- 
terials and the samples after synthesis was performed in 
an Ar filled glove box (M. Braun Ph 2 o < 0.1 ppm). 

It has been reported that a homogeneity range ex- 
ists for La2C3 extending from 45.2% to 60.2% atom-% 
C content^ After investigation of a series of La2C3_,5 
samples ; 13 ' 15 ' 16 it was concluded that samples with a car- 
bon deficit upon annealing phase separate into two su- 
perconducting phases with rather sharp T c 's of ~ 6 K and 
13.4 K. The latter transition which can be attributed to 
the almost stoichiometric La2C3 phase, is significantly 
higher than the T c of ~ 11 K reported so far for La2C3. 
Therefore, to obtain samples with a single sharp super- 
conducting transition at 13.4 K, it is essential to compen- 
sate possible losses of C in the arc melting and the sub- 
sequent annealing procedure by an excess of carbon up 
to 10%, although this may result in the impurity phase 
which is identified to be LaC2- Heat treatment of the 
sample buttons was performed in sealed Ta tubes un- 
der purified Ar atmosphere at 1000°C. After annealing 
at high temperatures, the samples were slowly cooled to 
room temperature with a rate of 5°C/hour. 

Neutron powder diffraction experiments on 2 samples 
of La2C3 (samples, SI and S2, cf. Table I) were per- 
formed at room temperature using the GEM difffactome- 
ter at the ISIS laboratory. For low temperature neu- 
tron powder experiments, another La2C3 sample (sam- 
ple S3) of ~ 10 g was sealed into a vanadium can un- 
der 1 bar of He gas and measured down to 5 K on the 
D2B diffractometer at the Institut Laue-Langevin (ILL). 
The structural parameters as well as the La2C3/LaC2 
composition were gained from two-phase Rietveld refine- 
ments using the FULLPROF package. 25 T c was deter- 
mined from dc magnetic susceptibility measurements us- 
ing a SQUID magnetometer (Quantum Design, MPMS 
XL magnetometer). The specific heat (C p ) was measured 
using a PPMS calorimeter (Quantum Design) employing 
the relaxation method. The C v contribution from the 
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LaC2 impurity phase was estimated from the results of 
a separate run on a pure LaC2 sample characterized also 
with neutron powder diffraction. To determine the upper 
critical fields, we measured the temperature dependence 
of the resistivity for a bar shape sample (lxl x 5 mm 2 ) 
of La2C3 at different magnetic fields up to 11 T. For 
magnetic fields higher than 11 T, we performed magneti- 
zation and resistivity measurements in a pulsed magnetic 
field up to 30 T at the Laboratoire National des Champs 
Magnetiques Pulses in Toulouse. 

III. NEUTRON POWDER DIFFRACTION 

La 2 C3 crystallizes with the cubic PU2C3 structure 
(I43d) with 8 formula units in the unit cell. As shown in 
Fig. [TJ C-C dumbbells are located in a distorted dodeca- 
hedral coordination ('bisphcnoid') formed by 8 La atoms. 
So far, several crystal structure determinations for La2C3 
have been carried out by neutron powder diffractio n 26 i 27 
and also by x-ray single crystal diffraction, 15 There was 
some variation found for the lattice parameters ranging 
from a = 8.804 to 8.818 A at room temperature^^^ 
Such a wide spread of structural parameters has been 
already known for Y2G^ as well as other C-containing 
superconductors like MgCNi3. 29 It is related to the C 
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FIG. 1: (Color online) Neutron powder diffraction patterns of 
La 2 C 3 (S3) at (a) 300 K and (b) 5 K (D2B, ILL). The (red) 
solid line represents the Rietveld refinement. The difference 
between the observed data and the calculated curve is shown 
at the bottoms while the Bragg positions are indicated by 
upper and lower bars for La2C3 and LaC2, respectively. The 
crystal structure of La2C3 and the C-C dumbbell structure 
surrounded by the La cage (bisphenoids) are shown in the in- 
sets of (a) and (b), respectively. Large (small) spheres denote 
La (C). 



content and often affects the superconducting properties 
quite substantially. For example, in MgC^Nis, the den- 
sity of state at Ep, N(Ep), decreases rather abruptly due 
to disorder-induced smearing of the electronic bands 1^ In 
addition, e-ph coupling mostly due to the low energy Ni- 
dominated phonon modes is reduced by C deficiency 31 ^ 32 
Similar effects of C deficiency may also cause the large 
variations of superconducting properties in -R2C3 com- 
pounds. 

Figure[T]shows the neutron powder diffraction patterns 
of La 2 C 3 (S3) collected at T = 300 K and 5 K. In addi- 
tion to the reflections which can be ascribed to La2C3, 
there are extra weak reflections, which belong to the im- 
purity phase, identified as LaC2 (I4/mmm). A two-phase 
Rietveld refinement was preformed to account for the ad- 
mixture of the LaC2 phase whose weight fraction is re- 
fined to ps 18%. The converged parameters of the La2C3 
phase include the lattice constants, the fractional coor- 
dinates (it, u, u)La of the La 16c site and (v, 0, l/4)c 
of the C 24c? sites, an isotropic (anisotropic) thermal pa- 
rameters for the La (C) sites. After refinements with the 
aforementioned parameters, in a last step, the C occu- 
pancy was varied. We checked convergence by varying 
the C occupancy in the cases where the other refined 
parameters are fixed or relaxed. For both C defi- 

ciency of ~ 2% is derived consistently. The results of the 
refinements for the two patterns collected at T = 300 K 
- 5 K are listed in Table U Also the results for samples 
SI and S2 obtained using time-of-flight (TOF) patterns 
collected at room temperature are given. 

The C-C distance at room temperature amounts to ~ 
1.295 A, in agreement with the previously reported value, 
1.296(9) A [Ref. 27] but smaller than that reported early 
on, 1.32 A. 2(5 Similar C-C distances have also been ob- 
served in binary rare earth dicarbides and ternary car- 
bide halides i 33 ! 34 ' 35 ! 36 From the studies of electronic be- 
havior on a series of the dicarbides, MC2, the C-C bond 
length inside the octahedral metal atom cage is found to 
be linked to the valence state of the M atoms^ In a sim- 
ple ionic picture, CaC2 (Ca 2+ (C2) 2 ~) with filled bonding 
7T states but empty antibonding ir* and a* states, has a 
short C-C bond distance of w 1.2 A, while with increasing 
electron count for the C2 unit e.g. in UC2 (U 4+ (C2) 4 ~), 
the C-C bond distance increases to ~ 1.35 A as the an- 
tibonding 7t* states are gradually filled. La2C3 has a 
shorter C-C bond in a C2 4_ unit, which is presumably 
due to the different degree of charge transfer to antibond- 
ing C-C 7r* states. These results indicate that the charge 
transfer between M states and antibonding tt* state of 
the C-C dimers strongly depends not only on the valence 
state of M, but also the local environment around the 
C-C dimers because of different cage structure. 

Interestingly we observed a significant temperature de- 
pendence of the C atom position in contrast to that of the 
metal atom cage which contracts in a regular way. With 
decreasing temperature a slight increase of the C-C bond 
distance is found (see Table UJ). A possible explanation 
of this observation could be a temperature dependence 
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TABLE I: Structural parameters for La2C3 obtained from Rietveld refinements of neutron powder diffraction patterns at T = 
300 K (Sf-S3) and at low temperature down to 5 K (S3), a is the cubic lattice parameter and dc-c is the C-C bond length, u 
and v are the fractional coordinates of the La 16c site (it, u, u)i,a and the C 24d sites (v, 0, l/4)c, respectively. The reduced x 2 , 
the occupancy for C are listed. The employed method, time-of-fiight (TOF) and constant-wavelength (CW) neutron powder 
diffraction is also indicated. 



Sample 


T (K) 


a (A) 


u (La) 


v (C) 


x 2 


dc-c (A) 


C occupancy 


Method 


SI 


300 


8.80991(8) 


0.05255(3) 


0.30155(9) 


2.074 


1.2942(16) 


0.983(8) 


TOF 


S2 


300 


8.8090(1) 


0.05268(6) 


0.30171(16) 


1.325 


1.2911(28) 


0.977(5) 


TOF 


S3 


300 


8.8096(3) 


0.05255(6) 


0.30144(14) 


8.96 


1.2977(19) 


0.982(2) 


CW 




200 


8.8013(3) 


0.05251(5) 


0.30114(14) 


9.42 


1.3001(17) 




CW 




100 


8.7936(3) 


0.05251(5) 


0.30085(13) 


10.06 


1.3041(17) 




CW 




5 


8.7904(3) 


0.05250(5) 


0.30065(13) 


10.29 


1.3071(17) 




CW 



of the charge transfer between the C-C ir* to the La d 
states. A carbon-metal orbital overlap will be affected 
by volume contraction at low temperatures, which will 
modulate the charge transfer between La d and C-C it* 
states. Another possible origin could be a tilting mo- 
tion of the C-C dimers leading to unusual non-elliptical 
temperature factors which can not be accounted for by 
the refinement procedure. Further studies are ongoing to 
clarify this issue. 

IV. ELECTRONIC STRUCTURE 
CALCULATIONS 

The electronic structure of La2C3 has been studied 
using the full-potential linear augmented plane wave 
(LAPW) method with local orbital extension, within the 
generalized gradient approximation (GGA) i 37 ' 38 i 39 We 
performed total-energy calculations for a series of the lat- 
tice parameters. For each lattice parameter the internal 
parameters (atomic positions) were relaxed according to 
the atomic forces. By this procedure, the fully relaxed 
structural parameters have been obtained including the 
lattice constant as well as the La and C atomic positions. 
The calculated equilibrium lattice constant is a = 8.829 
A, with a C-C distance of 1.333 A and a La-C distance 
of 2.709 A, corresponding to internal parameters, «La = 
0.0524 and vq = 0.2995. In comparison with the Rietveld 
refinement results at 5K, the calculated lattice parame- 
ter is within the normal error bar of LAPW calculations 
(~ 1%), while the internal atomic parameters for both 
La and C are in good agreement. 

Figure [2] displays the electronic structure with and 
without spin-orbit coupling included, using the optimized 
structural parameters. There are 6 C-C dimers in the 
primitive cell, and the orbitals on these form bonding and 
antibonding C-C 7r-7r* bands. The C 2s derived bands 
extend from -15.6 eV to -6.2 eV, relative to the Fermi 
energy. The lowest bands around -14 eV and the bands 
around -7.5 eV are the bonding and antibonding states, 
respectively. The bands around -4 eV are the bonding 
carbon p bands. The bands crossing the Fermi level, sep- 
arated by a 2.5 eV gap from the bonding C-C ir bands, 
are the hybridized La d and C-C antibonding n* states. 




L2 r h n r p 




12 r h n r p 



FIG. 2: Upper panel: Calculated band structure of La2C3 
with spin-orbit coupling, using the theoretical lattice constant 
with fully relaxed atomic positions. For comparison, the band 
structures near Ef with (middle panel) and without (bottom 
panel) spin-orbit coupling are also presented. 
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The electronic structure obtained including spin-orbit 
coupling is quite similar to that without spin-orbit cou- 
pling, except that the number of bands is doubled due 
to the asymmetric spin-orbit coupling effect (See middle 
and bottom panels in Fig. [2j. Lacking inversion sym- 
metry in the structure along with significant spin-orbit 
coupling, some spin-up and spin-down bands are mixed 
and the degeneracy is lifted. Therefore, in contrast to 
the previous conjecture^ the band splittings due to the 
asymmetric spin-orbit coupling do exist, and may also 
affect the superconducting properties. In particular, for 
La2C3, relatively flat bands are crossing Ep, thus such 
a band splitting may alter significantly the total elec- 
tronic density of states (DOS) at Ef. In Fig. [3J we show 
the DOS with and without spin-orbit coupling around 
Ep. The DOS near the Fermi level is characterized by 
a broad minimum centered at about 50 meV below E-p. 
As shown in the inset of Fig. [31 for stoichiometric com- 
pounds N(Ep) is reduced by ~ 10% by introducing the 
spin-orbit coupling. 

We also calculated the DOS for the system with 2% C 
deficiency by employing the virtual-crystal approxima- 
tion. We noticed a shift of peaks towards higher energies 
above Ep. As seen in Fig. [31 while the electronic struc- 
ture based on the stoichiometric sample puts Ep on the 
shoulder of the peak located near the Fermi level, Ep for 
the 2% C deficient sample lies close to a minimum. This 
leads to a reduction of N(Ep) by ~ 25% as compared 
to the DOS based on the stoichiometric La2C3 assuming 
no spin-orbit coupling. Introducing the spin-orbit cou- 
pling, the reduction of N(Ep) is even more pronounced 
and amounts to ~ 30 %. Thus slight C deficiency can 
significantly change the electronic properties. This pro- 
vides an explanation of the strong dependence of T c and 
the upper critical field on the composition found in Y2C3 
and LaaCgAia^i^ia 

Singh and Mazin concluded that X p h has a maximum 
at the stoichiometric band filling. 8 This is closely related 
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FIG. 3: (Color online) Electronic density of states of La2C3 
without (black) and with (yellow) spin-orbit coupling for 
the stoichiometric (solid) and 2% C-deficient (dashed) com- 
pounds. The inset shows the DOS near the Fermi level. 



to the observation that the Fermi energy falls on a peak in 
the electronic density of states. Hence, T c w 13.4 K in the 
present 2% C deficient La2C3 samples is already higher 
than the previously known T c ~ 11 K, but it seems that 
there is still a possibility to increase T c further for stoi- 
chiometric La2C3. In this view, the doping dependence 
of T c for La2C3 can also be understood. Th doping in 
La2C3 enhances T c up to 14.3 K^i while Lu or Y doping 
slightly decreases T c £2- Considering the electron count 
for Th, Th doping at the La sites will donate more elec- 
trons than La, in contrast to Y or Lu substitution. For 
the previously investigated La2C3 samples (T c ~ 11 K), 
the reduced electron concentration due to C deficiency, 
is compensated by Th doping, but not by Lu or Y dop- 
ing. Therefore, enhancement of T c by Th doping can be 
attributed to the recovery of the band filling towards the 
stoichiometric value. 



V. SPECIFIC HEAT 

Figure 2] shows the temperature dependence of the 
specific heat of La 2 C 3 (SI) at H = and H = 9 T. 
The contribution of LaC2 to the total sample capacity, 
which amounts to < 5% over the whole temperature 
range, was subtracted. Note that T c of LaC2 is < 1.6 
K , 42 ' 43 below the temperature range of the present mea- 
surements. There is no offset of C p /T at H= as the 
temperature approaches zero as can be seen in the in- 
set of Fig. [4] This proves that the C p contribution from 
the non-superconducting part, LaC2 has been completely 
accounted for the subtraction. A sharp anomaly at f=a 
13.4 K is clearly resolved indicating bulk superconduc- 
tivity in the La2C3 sample, which has not been observed 




100 150 
t (K 2 ) 



250 



FIG. 4: (Color online) Temperature dependence of the spe- 
cific heat of La2C3 at H= and 9 T. The solid (black) line 
through the data points at H = 9 T (also in the right inset) 
is a fit described in the text. The insets show the low temper- 
ature behavior of C P (T)/T and the magnetic susceptibility 
with zero-field cooled (ZFC) and field-cooled (FC) modes. 
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previously^ The onset of 13.4 K determined from the 
specific heat jump is consistent with that obtained from 
the susceptibility measurements. At H = 9 T, the super- 
conducting anomaly is shifted to lower temperatures. A 
clear deviation from a Debye T 3 law is observed in the 
normal state C' p , which can be attributed to a contribu- 
tion from low-lying Einstein phonon modest 

In order to account for the normal state C p , we fitted 
the data obtained at H = 9 T by a polynomial, 

C P (T) = 7N T + C latUce (T) = 7N T + [3T 3 + ST 5 . (1) 

Here, 7n is the Sommerfeld coefficient and (3 is related 
to the Debye temperature 0d (0) via 

03 (0) [K] = 1944n//3 [J/molK 4 ], (2) 

where n is the number of atoms per formula unit. Since 
the applied magnetic field H= 9 T is not sufficient to sup- 
press superconductivity completely, we used a constraint 
according to 

/ C p (T)/TdT= / (tn + f3T 2 + 5T 4 )dT, (3) 
Jo Jo 

to assure the entropy conservation for the superconduct- 
ing state. 

The solid line in Fig. H is the best fit to the H = 9 T 
data for T C (H = 9 T) < T < 15 K, yielding the param- 
eters, 7n = 10.60(4) mJ/mol K 2 , (3 = 228.2(11) /xJ/mol 
K 4 , and 5 = 0.5988 /iJ/mol K 6 . The measured 7n is 
much higher and more realistic than those previously 
reported from specific heat^ and upper critical fields 
measurements ^ The corresponding Debye temperature 




T(K) 

FIG. 5: (Color online) Superconducting part of the electronic 
specific heat. The (red) solid line is the best fit according 
to the a-model assuming an isotropic s-wave BCS gap as de- 
scribed in the text. The (blue) dashed line represents the 
BCS result. The inset shows the entropy conservation for the 
electronic specific heat in the superconducting state. 



©d(0) = 349(1) K is comparable with that found in pre- 
vious work.— 

The specific heat difference AC P (T)/T between the 
normal and superconducting state and the entropy for 
the superconducting state is shown in Fig. [5] For com- 
parison, we also plot the BCS curve for the weak e-ph 
coupling limit. A clear deviation from the BCS curve is 
observed. The C p anomaly at T c as well as the inter- 
section temperature where AC P /T = is much higher 
than expected from the BCS prediction. The solid (red) 
line is the theoretical fit based on the 'a-model' assum- 
ing an isotropic s-wave BCS gap A(T) scaled by the ad- 
justable parameter, a — A(0)//cbT c 4£ For the weak cou- 
pling limit a is 1.76. The detailed temperature depen- 
dence of ACp/T was fitted by two adjustable parameters: 
a and the Sommerfeld coefficient y^. The data are very 
well reproduced by a = 2.43 and 7n = 8.97 mJ/mol K 2 . 
The 7n value from the a model fit is slightly smaller than 
obtained from the normal state C p , but in view of the 
experimental resolution and the accuracy of the model, 
it can be considered as satisfactory agreement. The nor- 
malized specific heat jump, AC p /^T c is also higher than 
the weak limit BCS value, 1.426 giving a clear indication 
for an enhanced e-ph coupling. 

The normalized electronic specific heat, C es /^T C in 
the superconducting state is shown in Fig. [S] The solid 
line is an exponential fit to the data for 2.0 < T c /T < 
4.5 using the form C es /^T C oc exp(-0.82aT c /T) with 
a =2.31(5). C es exponentially vanishes for T — ► K, 
clearly manifesting the absence of gap nodes in the su- 
perconducting order parameter. The a value is somewhat 
lower than found from the a-model fit in which the ra- 
tio A(0)//cbT c is largely determined by the shape of the 




TIT 

c 



FIG. 6: (Color online) The electronic contribution of the spe- 
cific heat Ces is plotted on a logarithmic scale versus T c /T. 
The (red) straight line is an exponential fit for 2.0 < T c /T < 
4.5. The magnetic field dependence of Sommerfeld coefficient 
y(H) is presented in the inset with the dashed line as a guide 
to the eye. 
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C p jump near T c . However the discrepancy is less than 
10%, which is probably due to a restricted temperature 
range. The magnetic field dependence of j(H) at T — > 
K is estimated from the linear fit of C p /T versus T 2 
for 2 K < T < 4 K, and it is found to increase linearly 
with H for low fields up to H ss 0.4 H C 2. This behavior 
consistently supports a fully-gapped and almost isotropic 
superconducting order parameter 4i 

The thermodynamic critical field H C (T) can be deter- 
mined using H C (T) = y/SnAF where AF is the free 
energy extracted from the specific heat in the supercon- 
ducting state, AC P = -Td 2 (AF)/dT 2 . The temperature 
dependence of H C (T) is shown in Fig. [Jj Based on the 
BCS theory, H c (0) is derived from a fit to the tempera- 
ture dependence of H C (T) for T — > which is given by 



Y2C3. 48 Since the Fermi surface consists of several sheets 
from hybridized La d and C p states, we cannot rule out 
such a possibility for La2C3 as well. However, as shown 
above, the specific heat for the superconducting state can 
be well explained in terms of a single superconducting 
gap. Usually, in a system where the d bands are hy- 
bridized with s or p bands, the disparity between the 
bands is reduced, and impurity scattering significantly 
smears out the two-gap superconductivity ! 49 ' 50 Consid- 
ering the polycrystalline nature of the sample, therefore, 
it is rather unlikely that, if any, two-gap superconductiv- 
ity survives. Further studies are required to clarify the 
possible two-gap superconductivity in the sesquicarbidc 
systems. 



Hc(T) 
H c {0) 



= 1 



- - ( y 



(4) 



where the empirical parameter (3 is ~ 1 for the weak 
coupling limit but reduced with increasing e-ph cou- 
pling. The solid line in Fig. [7] is a best fit yielding 
H c (0)= 151.3(1) mT and /?= 0.7933(3) indicating as 
well strong e-ph coupling. The deviation function, D(t) 
= H c (T)/H c (0)-(l-t 2 ) with t = T/T c as well provides 
clear signatures for strong e-ph coupling. In the case of 
strongly e-ph coupled superconductors, D(t) is positive 
and shows a broad maximum while the weak-coupling 
BCS curve exhibits a negative dip at t 2 « 0.5. As shown 
in the inset of Fig. [JJ D(t) for La 2 C3 passes through a 
maximum at t 2 m 0.5. 

Based on the specific heat results for the superconduct- 
ing states, we conclude that La 2 C 3 is an s-wave, single 
gap superconductor with strong e-ph coupling. Recently, 
Harada et al. from 13 C nuclear- magnetic-resonance mea- 
surements suggested a multi-gap superconductivity for 




I" 60 -X 



T(K) 

FIG. 7: (Color online) The temperature dependence of the 
thermodynamic critical field H C (T) estimated from the spe- 
cific heat in the superconducting state. The solid (red) line 
is a best fit (see the text). The inset shows the deviation 
function, D(t) = H C (T)/ H c (0)-(l-t 2 ) with t 2 = (T/T c f. The 
weak coupling BCS curve is also presented for comparison. 



VI. UPPER CRITICAL FIELDS 

Figure [8^a) shows the temperature dependence of the 
resistivity under different magnetic fields up to H = 11 T. 
T c is determined at 50% decrease from the normal state 
resistivity value and the transition width is taken as the 
temperature interval between 10% and 90% of the transi- 
tion. In order to determine H&iT) at low temperatures, 




FIG. 8: (Color online) (a) Temperature dependence of the re- 
sistivity (p(T)) under magnetic fields from to 11 T increas- 
ing by AH = 1 T, from right to left, (b) Magnetoresistivity 
(p(H)) at different temperatures at T = 10,..., 2.4 K from left 
to right. The (black) solid line presents 50% of the resistive 
transition, while two (black) dashed lines above and below de- 
note 90% and 10% of the resistive transition. Note that p(H) 
in (b) shows a slight positive magnetoresistance. The inset of 
(b) presents the first derivative of the magnetization versus 
field curve at T = 4.2, 2.5, and 1.4 K, from left to right. 
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the magnetoresistance was measured up to w 30 T at 
different temperatures from T = 1.8 to 10 K (Fig.^b)). 
The transition width remains almost unchanged down to 
low temperatures. As a consistency check we also carried 
out high field magnetization measurements at low tem- 
peratures. H C 2(T) is determined as a kink of the first 
derivative of M (H) (See the inset of Fig. |5J|. 

The H C 2(T) data obtained from the three different 
methods are compiled in Fig. [9] For comparison, we plot- 
ted the Werthammer-Helfand-Hohenberg (WHH) pre- 
diction for conventional superconductors. 51 There is a 
clear deviation from the WHH behavior in the H C 2(T) 
curve. This is in contrast to the previous report on 
the H c2 {T) of La 2 C 3 sample (T c ~ 11 K), which fol- 
lowed the WHH prediction rather well. Such an enhance- 
ment of H C 2(T) has been attributed to several different 
origins such as localization effects in highly disordered 
superconductors^ anisotropy of the Fermi surface^ 
and strong e-ph coupling ] 54 ' 55 First, considering disor- 
der effects, the application of the magnetic field weakens 
the localization effects and thus reduces the Coulomb 
pscudopotential (/J.*). This in turn strengthens the su- 
perconductivity and leads to an enhancement of H C 2 as 
the temperature decreases and larger fields have to be 
applied to drive the system into the normal state. There- 
fore, the enhancement of H C 2 is closely linked to the neg- 
ative magnetoresistance in the normal stated The pos- 
itive magneto-resistance observed in La2C3 (See Fig. [8]), 
however, indicates that the localization effects tend to 
be enhanced under high magnetic fields, which rules out 
the possibility of disorder- induced enhancement of H C 2- 
Secondly, when the FS is distorted from the spherical 
shape, an increase of H C 2 becomes more pronounced due 
to an anisotropy in the FSi 5 ^ Electron structure calcu- 
lations for La2C3 unveil rather complex multisheets of 
the Fermi surface, similar to that of Y2C3 (cf. Fig. 5 
in Ref. H). A quantitative comparison between experi- 
ment and theory is beyond the scope of this work, but 
considering the distorted FS shape, the FS anisotropy 
will most likely contribute to the enhancement of H C 2 at 
low temperatures. Note, however, that the previous re- 
port on H c2 {T) for La 2 C 3 (T c ~ 11 K) with C deficiency 
showed good agreement with the WHH curved With C 
deficiency, the detailed FS can be modified, but the over- 
all FS shape would be preserved, thus similar enhance- 
ment of H C 2 at low temperatures is expected even for C- 
deficient samples with lower T c . Therefore it seems that 
the FS anisotropy cannot be the only source for H c2 en- 
hancement. Another possibility is strong e-ph coupling. 
When the e-ph coupling is in the strong coupling regime 
(Xph > 1), H C 2(T) starts to deviate from the WHH curve 
and remains linear down to lower temperatures showing 
finally an upward curvature for higher e-ph coupling. 55 
As shown in Fig. [51 H C 2{T) for La2C3 shows an almost 
linear temperature dependence down to ~ 0.2T C . From 
the specific heat studies, we already found evidence that 
the e-ph coupling is in the strong coupling regime, thus 
it will also affect the temperature dependence of H c2 for 



La 2 C 3 . 

From an extrapolation, we obtained H C 2(0) ~ 19 T. 
Even though H C 2(0) is much more enhanced than the 
WHH prediction, it is still clearly below the Pauli limit 
of H p (0) = A(0)/2 v / 2^b = 1-83 k B T c w 24.5 T for the 
weak coupling limit. Considering the strong e-ph cou- 
pling with enhanced A(0)/k B T c = 2.31 - 2.43, H p (0) can 
be even higher up to ~ 32 T. Therefore, H c2 (0) w 19 T, 
well below the Pauli limit, indicates that H C 2 is mainly 
determined by orbital depairing under magnetic fields. 
This result suggests that the possible effect due to non- 
centrosymmctry is not significant for H C 2 as we will dis- 
cuss below (see Sec. VII). 

Based on the H C 2(0) and H c (0) values, we estimated 
the superconducting parameters for La2C3. According 
to the Ginzburg-Landau (GL) theory, the upper critical 
field H C 2 (0) , the lower critical field H c i (0) and the ther- 
modynamic critical field H c (0) can be described with the 
GL coherence length £(0) and the GL parameter re(0) = 
A(0)/f(0) according to 



H c2 (0) 



2<(0)' 



H cl (0) 



g c 2 (Q) 

H c2 (0) 



[ln/c(0) + 0.08], 



H c (0) = 



H c2 (0) 
\/2k(0)' 



(5) 



(6) 



(7) 



where $0 is the flux quantum. With H C 2(0) = 19 T and 
H c (0) = 151 mT, we obtained £(0) « 42 A and A(0)= 373 




2 4 6 8 10 12 14 
7(K) 



FIG. 9: (Color online) The temperature dependence of the 
upper critical field for La2C3 samples, SI (square) and S3 
(circle), estimated from p(H) (filled symbols), M(H) (crossed 
symbols), and p(T) (open symbols) measurements. For com- 
parison, the Werthammer-Helfand-Hohenberg curve with the 
Maki parameter a = and the spin-orbit scattering Aso = 
(see Ref. l5lh is presented with (blue) dashed line. The Pauli 
limit H p is indicated by the arrow. 
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nm. The GL parameter, k(0) is w 90, indicating type- 
II superconductivity. Similar values have also been re- 
ported for the 14 K phase of Y 2 C 3 (T c = 13.9 K): H cl (0) 
= 3.5 mT, A(0) = 430 nm, H c2 = 24.7 T, and £(0) = 36 



VII. DISCUSSION 



rithmic averaged phonon frequency: 
2A(0) 



3.53 



AC P (T C ) 



1 + U.bx'^ln- 



= 1.43 



1 + 53aT 2 ln~ 



(8) 



(9) 



First we briefly discuss the effect of non- 
centrosymmetry in La2C3. There have been several 
characteristic features reported concerning the effect 
of the lack of centrosymmetry on the superconducting 
properties : (1) enhancement of H C 2(0) above the Pauli 
limit (2) presence of a line node in the superconducting 
gap parameters. Furthermore, it has been suggested that 
a description of the superconducting properties should at 
least include two bands due to the non- centrosymmetry 
in rare earth sesquicarbidesi 2 ^ Our electronic structure 
calculations confirm the splitting in the electronic bands 
near E-p indicating that the spin degeneracy is lifted due 
to the sizable asymmetric spin-orbit coupling. 

However, even though H C 2(0) is clearly enhanced, it 
does not exceed the paramagnetic limit (Fig. [9]). The 
specific heat at low temperatures also reveals that the su- 
perconducting gap in La2C3 has an isotropic s-wave sym- 
metry. Therefore the effect of the non-centrosymmetry 
appears to be not significant in La2C3. Note that the 
band splitting with respect to the superconducting tran- 
sition temperature is much smaller for La2C3 than in 
other non-centrosymmetric superconductors (See Table 
ITT)) . In Li2Pd3B where the band splitting due to non- 
centrosymmetry is comparable with that of La2C3, con- 
ventional BCS type behavior with the isotropic super- 
conducting gap has been found from penetration depth 
measurements using muon-spin-rotation experiments. 59 
Therefore for La2C3, the asymmetric spin-orbit coupling 
appears to be not strong enough to induce a significant 
effect on the superconducting properties. In this respect, 
it would be very interesting to study heavier rare earth 
metal carbides such as LU2C3 (T c w 15 K)££ and TI12C3 
(T c » 4 K)4i 

For e-ph coupled superconductors, CarbotteSi has pro- 
posed that the characteristic thermodynamic quantities 
follow empirical formulas which can be described by one 
adjustable parameter, x — u>i n /T c where w;„ is the loga- 



TABLE II: Energy splitting of the bands at E F (AE) and T c 
for several non-centrosymmetric superconductors. 





AE (meV) 


AE/k B T c 


Tc (K) 


Ref. 


La2C3 


20 - 30 


~ 20 


13.4 


this work 


CePt 3 Si 


50 - 200 


> 1000 


0.75 


21 


Cd 2 Re 2 07 


~ 70 


~ 700 


~ 1 


57 


Li 2 Pt 3 B 


~ 200 


> 1000 


~ 2 


58 


Li 2 Pd 3 B 


20 - 50 


> 35 - 85 


~ 7 


58 



AC P (T) - AC P (T C 
JnT c - 7 N T 



-3.77 



1 + 117a;- 2 ln 



2.9 



7nT c 2 

#e 2 (0) 



0.168 



1 - 12.2ar 2 ln- 



H c (0) 



dH c (T)/dT\ T T c 



= 0.576 



1 - 13.4.x- 2 ln— 
3.5 



(10) 



(11) 



(12) 



This analysis has been successfully applied to various 
metal-alloy superconductors^ and also to other re- 
cently discovered carbon-contained superconductors such 
as borocarbides£2i22 and MgCNi 3 ^i 

Figure [TO] shows the thermodynamic quantities with 
variation of x — uJi n /T c according to Eqs. (8)-(12). The 



5.5 
5.0 
4.5 
4.0 
3.5 
2.8 

2.4 
2.0 

1.6 
-4 
-8 
-12 
-16 
-20 
0.16 



2A(0)//cT 




(AC(7)-AC(T))/( TN r-Y N T) 
H 1 1 1 1 ' 1 ' 




FIG. 10: (Color online) Several thermodynamic quantities 
as a function of the logarithmic averaged phonon frequency 
(solid line) according to Eqs. (8)-(12). The thermodynamic 
quantities estimated from the specific heat data for the sam- 
ples are indicated with (red) open circles for SI, (black) 
crosses for S2 and (blue) filled squares for S3. 
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measured values extracted from the C p data are also 
plotted onto the empirical curves. Each thermodynamic 
quantity provide a value for x — ui n /T c as indicated for 
three La2C3 samples (S1-S3) in Fig. [TO] From five x's 
and T c = 13.4 K, we obtain a mean value uj[ n = 109 K 
± 24 K, 101 K ± 24 K, and 121 K ± 32 K for the sam- 
ple SI, S2, and S3, respectively. Averaging all w;„'s for 
the samples we conclude on ioi n — 110 K ± 27 K. This 
uji n value is only ~ 30% of the Debye temperature, wd 
w 350 K. The Debye frequency is determined only from 
the phonon DOS profile, F(ui), while LO\ n results from a 
weighting by the e-ph coupling function a 2 {uj). Such a 
reduced u>i n value compared to cjd indicates the impor- 
tance of the low energy phonon modes for the supercon- 
ductivity. The deviation of the normal state C p from the 
Debye T 3 law (See. Fig 2]) supports the presence of the 
low-lying Einstein phonon modes. 

From ab-initio calculations for Y2C3, the frequencies 
of the full symmetry Raman modes are estimated to be 
175 cm -1 for a Y dominated mode and 1442 cm -1 for an 
almost pure C-C bond stretching mode.— The relevant 
phonon modes for the superconductivity, therefore, are 
the La dominated phonon modes rather than the high 
frequency C-C bond stretching modes. From simple con- 
version based on the atomic mass difference between La 
and Y, we can estimate the frequency of the La domi- 
nated Raman mode to be ~ 140 cm" 1 . In comparison 
with LOi n m 76 cm -1 ±15 cm" 1 , this phonon frequency is 
somewhat higher. However, since the lattice constant of 
La2C3 is larger than that of Y2C3, the interatomic forces 
are expected to be somewhat weaker for La2C3. The es- 
timated frequency of the La dominated phonon mode, 
therefore, is rather overestimated, and it can be reduced 



30 
25 
20 



* 15 



10 
5 



La 2 c 3 


/ D =349(1 ) K 
'' u* = 0.1 

^^^^^^^^^^^^ 




«a > =110 ±26 K 

In 

u* = 0.1 



0.5 



1.0 



1.5 



ph 



2.0 
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FIG. 11: The calculated T c from the modified McMillan for- 
mula (solid line) as a function of the e-ph coupling constant, 
X p h- The borders of the shaded (grey) areas represent the 
results obtained for the upper and lower bounds of uji n . For 
comparison, we also plot the calculated T c from the origi- 
nal McMillan formula with dashed line. The horizontal (dark 
grey) solid line denotes the measured T c . 



and come closer to the u>i n . Recent ab-initio calculations 
for Y2C3 also reported that most of the e-ph coupling 
emerges from the low frequency Y related phonon modes, 
not from the high frequency C-C phonon modest 

Other relevant phonon modes for superconductivity in 
La2C3 could be the tilting vibrations of the C-C dimer. 
For AC2 (A = Ca, Sr, and Ba) containing the C-C dimers 
inside the octahedral metal atom cage, it has been shown 
that the binding of the C-C dimer to the surrounding oc- 
tahedral metal atom cage is weak, thus allowing a struc- 
tural transition from tetragonal to a cubic phase at high 
temperatures as well as several modifications of the C-C 
dimer structure at low temperatures. 66 Raman scattering 
studies on rare-earth carbide halides also showed that the 
frequency of the C-C tilting modes is ~ 400 cm" 1 , much 
lower than the stretching modes with a frequency of ~ 
1590 cm" 1 .— Compared to the octahedral metal atom 
cages, the size of a bisphenoid La cage is larger, thus 
the binding of the C-C dimer to the La cage could be 
even weaker. Consequently the tilting modes of the C-C 
dimer in La2C3 are expected to have lower frequencies, 
which can allow sizable contribution to e-ph coupling. 
Further studies on the C-C dimer tilting modes e.g. Ra- 
man scattering experiments or ab-initio calculations are 
highly desirable. 

Finally we estimate the e-ph coupling constant for 
La2C3. Based on the Sommerfeld coefficient and the 
electronic structure calculation, we obtain X p h from 
the determination of 7n using the equation 7n = 
{2n 2 kl/3)N(E F )(l+\ ph ). From the DOS calculations 
for the 2% C deficient compounds, we take N(Ep)= 1.75 
states/eV f.u. Using 7n =10.6 mJ/molK 2 , the estimated 
X p h is 1.4 indicating that the La2C3 clearly belongs to the 
strong coupling regime. 

In Fig. [11] we plot the calculated T c from the modified 
McMillan formula, 



L2 



exp 



-1.04(1 + Ap h ) 



\ ph - (1 + 0.62X ph )n* 



(13) 



with the measured T c = 13.4 K (grey horizontal line). 
For comparison we also plot the calculated T c from the 
original McMillan formula^ In both cases, the Coulomb 
pseudopotential fi* is fixed to 0.1. With the original 
McMillan formula assuming the effective phonon fre- 
quency to be the same as the Debye frequency, we found 
that the measured T c is reproduced with weak e-ph cou- 
pling, Xph ~ 0.75. However such a weak coupling cannot 
reconcile with the results of specific heat and upper crit- 
ical fields as discussed above. If we use the logarithmic 
averaged phonon frequency toi n in the modified McMillan 
formula, Eq. (13), a large value of X p h is necessary to re- 
produce the measured T c . Because of the relatively large 
error for the calculated T c from the modified McMillan 
formula due to the reduced accuracy for determining wi n , 
it is difficult to determine the X p h value precisely. How- 
ever it is obvious that X p h is larger than 1.3, which is in 
good agreement with the estimate from the comparison 
of the Sommerfeld coefficient. 
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In conclusion, specific heat and upper critical field 
studies were performed on the sesquicarbide supercon- 
ductor La2C3 together with neutron powder diffraction 
as well as electronic structure calculations. The main 
conclusions are the followings. (i) The density of states 
near the E-p is very sensitive to a small C deficiency 
in the sample, which provides an explanation for the 
wide scatter of T c 's and superconducting properties ob- 
served in the previous reports, (ii) The temperature and 
magnetic field dependence of the specific heat is consis- 
tent with a single gap s-wave BCS superconductor, (iii) 
The logarithmic averaged phonon frequency for the su- 
perconductivity is quite low suggesting the importance 
of low energy phonon modes for the superconductivity, 
(iv) Even though the band splitting due to the non- 
centrosymmetry in the structure and the spin-orbit cou- 
pling is clearly confirmed by the electronic structure cal- 



culations, its effect seems not sizable enough to cause ex- 
otic superconducting properties as observed in other non- 
centrosymmetric superconductors. To conclude, all of the 
summarized features suggest that La2C3 is a strongly 
coupled BCS-type superconductor with an isotropic s- 
wave superconducting gap. 
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